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Abstract. High- resolution spectroscopic time series of rapidly oscillating Ap stars show 
evidence for a co-existence of standing and running waves in their atmospheric layers. With 
the purpose of understanding these observations we have carried out a theoretical analysis 
of the pulsations in the outermost layers of these stars, starting from the simplest possible 
model that still retains all important physical ingredients. In our analysis we considered an 
isothermal atmosphere in a plane-parallel approximation. Moreover we assumed that in the 
region considered the magnetic pressure is much larger than the gas pressure and, consequently, 
that the magnetoacoustic wave has decoupled into its acoustic and magnetic components. Using 
the analytical solutions for the velocity components appropriate to this model we estimate the 
velocity component parallel to the line of sight averaged over the visible stellar disk. Fitting the 
latter to a function of the form Acos(crt+phase), with a the dimensionless oscillation frequency 
and t the dimensionless time, we derive the amplitude A and the phase for our model as function 
of height in the atmosphere. 



1. Introduction 

The rapidly oscillating Ap (roAp) stars are cool chemical peculiar stars which are located in 
the main sequence part of the classical instability strip. They pulsate with periods that vary 
typically from 5 to 21 minutes, e.g. [1], and have strong large scale magnetic fields, with typical 
intensities of a few kG. Moreover they have masses of M=2M Q [2]. Due to these characteristics, 
the roAp stars offer the opportunity to observe the interaction of p-modes with strong magnetic 
fields as can be done for no other star but the sun. The scientific interest of these stars has been 
recognised since the discover of their oscillations and, in the past few years, a several number of 
observational results have been published. High-resolution spectroscopic studies done for these 
stars show evidence for a co-existence of standing and running waves in their atmospheric layers, 
e.g. [1], [3], [4], [5]. 



2. Theoretical Modeling 

With the purpose of understanding these observations we have carried out a theoretical analysis 
of the pulsations in the outermost layers of these stars, starting from the simplest possible 
model that still retains all important physical ingredients. In our analysis we considered an 
isothermal atmosphere in a plane-parallel approximation. The photosphere of a roAp star can 
be located either in the magnetoacoustic, or in the magnetically dominated region, depending on 
the strength of the magnetic field, and, due to that, depending on the region that is considered, 



the oscillations in the atmosphere of these stars might look significantly different [6]. In this 
study we assumed that the region considered is magnetically dominated, i.e., the magnetic 
pressure is much larger than the gas pressure and, consequently, that the magnetoacoustic wave 
has decoupled into its acoustic and magnetic components. The magnetic field is considered to 
have a dipolar geometry and, hence, to vary in intensity and inclination with latitude. All the 
results presented are for a model of mass M=2.0 M Q and radius R=2.0 ~Rq, Mq and being, 
respectively, the mass and the radius known for the sun. The temperature of the isothermal 
atmosphere is T=8440 K and the atmosphere is matched continuously onto the polytropic 
interior at a pressure P=7563 gcm _1 s~ 2 , a value motivated by a stellar model with a similar 
mass and radius. The cyclic oscillation frequency v is assumed to be 3.08 mHz, which is above 
the critical frequency in this atmosphere, and the characteristic magnetic field is assumed to 
have a strength of 3000 G. The observer is considered to be pole-on in relation to the magnetic 

field axis. In Fig.l we show a schematic representation of the star, where B is the local magnetic 
field, a is the angle between the local magnetic field and the direction of the observer, 9 is the 
co-latitude and <f) the longitude. 



Starting from the equations for the displacement parallel and perpendicular to the direction 
of the magnetic field in the region where the magnetic pressure dominates over the gas pressure 
[7] and using the analytical solutions for the velocity components appropriate to this model we 
determine the expression for the velocity component parallel to the line of sight averaged over 
the visible stellar disk, V{ n t, 
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Fig.l - Schematic representation of the star. B is 
the local magnetic field, a is the angle between 
the local magnetic field and the direction of 
the observer, 6 is the co-latitude and cj> the 
longitude. 
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and A c (9) is the amplitude of the acoustic wave, a is the dimensionless angular frequency, 
a the limb darkening coefficient, p the dimenssionless pressure, (j) c (0) the phase of the acoustic 
wave, 77 the dimenssinless depth, A m (9) the amplitude of the magnetic wave, 7 the first adiabatic 
exponent, b the magnetic field modulus, H the dimensinless scale height, p the dimensionless 
density, (j) m (6) the phase of the magnetic wave and C s a constante that is related to the sound 
speed in the atmosphere. 

Fitting Vi n t to a function of the form Acos(crt+</>), we derive the phase (f>, 
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and the amplitude A, 




for our model as function of height in the atmosphere, where C c is defined as 
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3. Results 

In Fig. 2 we show the amplitude and phase diagrams as function of height in the atmosphere. 
The results presented assumed an A c (8) and an A m {9) varying according to, 

A c {9) = 0.1 cos 9, 



A m {9) = sine, 

These functional form for the amplitudes were motivated by the results Sousa and Cunha [7]. 
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Fig. 2 - Amplitude and phase diagrams as function of heigth in atmosphere, z is the vertical 
local coordenate, defined to be zero at the photosphere and increase outwardly.In the left panel 
we show the amplitude and phase derived when the velocity is integrated over the whole stellar 
disk, < 9 < 7r/2, and in the right panel the amplitude and phase derived when the velocity is 
integrated over 7r/4 < 9 < tt/2. 



In the left panel we show the amplitude and phase derived when the velocity is integrated 
over the whole stellar disk, < 9 < tt/2, and in the right panel the amplitude and phase derived 
when the velocity is integrated over 7r/4 < 9 < tt/2. The second case is considered due to the 
fact that the observed spectral lines, from which the velocity is determined, might originate in 
particular regions of the stellar surface, rather than being homogenously distributed over the 
whole surface. It can be seen that in both cases the amplitude shows an overall increase in the 



outward direction. This is the result of the decreasing density with height. At first sight, the 
phases show a pattern that seems to indicate the existence of an outwardly running wave, when 
the phase is decreasing with height, an inwardly running wave, when the phase is increasing 
with height, and a standing wave when the phase is almost constant. But, in fact, all of these 
patterns are the result of summing two waves, one acoustic outwardly running wave and one 
magnetic standing wave, in different proportions. The amplitudes and phases derived from 
v^nt, have contributions from the acoustic and magnetic components and depending on which 
is dominant different results are found. Thus, when the acoustic contribution is dominant we 
see a running wave propagating outwardly, when the magnetic contribution is dominant we see 
a standing wave, and when the two waves have similar contributions, the phase is such that 
it might look like a wave propagating inwardly. So, with this study we suggest that what is 
observed in high-resolution spectroscopic time series of roAp stars results from the superposition 
of different contributions of acoustic and magnetic waves. 
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